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Abstr act. We introduce a family of weight matrices W of the form T(t)T* (t), 
T{t) = e^*e®* , where is certain nilpotent matrix and 3 is a diagonal 
matrix with negative real entries. The weight matrices W have arbitrary size 
N X N and depend on A'^ parameters. 

The orthogonal polynomials with respect to this family of weight matrices 
satisfy a second order differential equation with differential coefficients that are 
matrix polynomials F2, Fi and Fo (independent of n) of degrees not bigger 
than 2, 1 and respectively. 

For size 2x2, we find an explicit expression for a sequence of orthonor- 
mal polynomials with respect to W. In particular, we show that one of the 
recurrence coefficients for this sequence of orthonormal polynomials does not 
asymptotically behave as a scalar multiple of the identity, as it happens in the 
examples studied up to now in the literature. 



1. Introduction 

In the last few years a large class of families oi N x N weight matrices W having 
symmetric second order differential operators of the form 



(1.1) [j^) ^^w+UJ 

has bee n introd uced (see [DGT] . [DG3] . [D2] . [D3] . [DdiT] . [DdI2) . [G], [GPTlj . 
|GPT2) . |CMV| . [Ft] ). The differential coefficients F2, Fi and Fq are matrix poly- 
nomials (which do not depend on n) of degrees less than or equal to 2, 1 and 0, 
respectively. As usual, the symmetry of an operator D with respect to the weight 
matrix W is defined by / D{P)dWQ* — J PdW{D{Q))* , for any matrix polyno- 
mials P, Q. 

A sequence {Pn)n of orthogonal polynomials with respect to a weight matrix W is 
a sequence of matrix polynomials satisfying that P-m n > 0, is a matrix polynomial 
of degree n with non singular leading coefficient, and J PndWP*^ — AnSn,m, where 
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A„, rt > 0, is a positive definite matrix. When A„ = /, we say that the polynomials 
{Pii)n are orthonormal (we denote by / the identity matrix). 

Just as in the scalar case, any sequence (P„)„ of orthonormal polynomials with 
respect to a weight matrix satisfies a three term recurrence relation 

(1.2) tP„(<) =A,+iP„+i(f)+B„P„(t) + A;P„_i(i), n>0, 
where v4„, n>l, are nonsingular and P„, ?i > 0, Hermitian. 

If (Pn)n is a sequence of orthonormal matrix polynomials with respect to W ^ the 
symmetry of the second order differential operator (jl.ip is equivalent to the second 
order differential equation 

(1.3) Pn{t)F2{t) + P,;(i)Pi(t) + Pnit)Foit) = A„P„(t), 
where A„ are Hermitian matrices (see Lemma 4 of |D1| ). 

The theory of matrix valued orthogonal polynomials was started by M. G. Krein in 
1949 [K1[|K2] (see also [Be) or |Atk| ). but more than 50 years have been necessary 
to see the first examples of orthogonal matrix polynomials satisfying that kind of 
differential equations (see |DG1[ [Gl IGPTlj ) . These examples will likely play, in the 
case of matrix orthogonality, the role of the classical families of Hermite, Laguerre 
and Jacobi in the case of scalar orthogonality. 

As their scalar relatives, these families of orthogonal matrix polynomials also satisfy 
many formal properties, relationships and structural formulas (see |DG2i[DLllD4| ). 
These structural formulas have been very useful to compute explicitly the orthonor- 
mal polynomials related to several of these families, and, in particular, their recu- 
rrence coefficients (|1.2p . The recurrence coefficients of these examples asymptoti- 
cally behave as scalar multiples of the identity. More precisely, either lim„ = al 
and lim„ P„ = bl f |DLl IGdI| ). or lim„ a„A„ = a/, a 7^ 0, and lim„ 5„P„ = 6/, for 
certain divergent sequences (a„)„, (6„)„ of positive real numbers f [DG21 iDLl ID4) ). 

The purpose of this paper is to introduce a new family of weight matrices ha- 
ving orthonormal polynomials satisfying second order differential equations and 
whose recurrence coefficients do not asymptotically behave as scalar multiples of 
the identity. 

Our example is the family of weight matrices of arbitrary size N x N constructed 
from the — 1 non-null complex parameters ai, ■ ■ ■ a^^i and the positive real 
parameter 6 (6 7^ 1) as follows. Consider the nilpotent matrix A and the diagonal 
matrices ^ and 5* defined by 

\ /O ••• \ 

1 ••• 

: ^ ^ 2 ••• 

aN-i : : : ■. : 

0/ \0 00---A^-l/ 



(1.4) 



(1.5) 



A 



^ = I 



/O ai 



\0 

6-1 





02 






N 



l^- 
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Let the diagonal matrix ^ and the upper triangular nilpotent matrix si be defined 

by 

(1.6) ^ = _^vi/-i, 

[f]-l 

(1.7) E <^,A^'^\ 

j=o 

where 

_ (l-b)^(2j + l)J--i 
^" (46)^(7V-l)Jj! ' 
Our weight matrix W is then defined by 

(1.8) W(t) = T{t)T*{t), T{t) = e^*e®*'. 

Since jz/ is nilpotent of order N ^ e^* is always a polynomial of degree iV — 1. 

For 6=1 (considering = 1), we recover the example 5.1 of |DG1) . 

For the benefit of the reader, we display here our weight matrix for size 2x2. For 
= 2 we have 



a \ ^ / -I 



and then 



(1.9) T= -^'^^ 

y ate"* e'* J \ ^ ^ J 

The content of this paper is the following. In Section 3, we prove that our weight 
matrix W always has a symmetric second order differential operator like (II. ip : 

Theorem 1.1. The second order differential operator U.l\) with differential coeffi- 
cients F2, Fi and Fq given by 

(1.10) ^^2(t)-* + ^p^K^]t, 

(1.11) Flit) = 2s/-i' + 2 (-bl + ^]) t, 



N - 1 

2, 



(1.12) Foit) ^2bJ 

is symmetric with respect to the weight matrix W il.8\) ( as usual [X, Y] denotes the 
commutator of the matrices X, Y : [X, Y] = XY — YX ). 



For N = 2, these differential coefficients are 



The rest of the Sections are devoted to study in depth the orthogonal polynomials 
with respect to our weigh matrix for size 2x2. The study of the orthogonal 
polynomials for higher size iV, iV > 3, remains a challenge. 
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In Section 4, we construct the following Rodrigues' formula for a sequence of or- 
thogonal polynomials with respect to the weight matrix W (II. 9|) : 



Theorem 1.2. Let the function Pn, n > 1, be defined by 
(1.13) 

. lap ^ - 



2t + e* y/^n[ Erf(V&i) - Erf(t) 



at 
2 



where Erf denotes the error function Erf(z) ~ —= e ^ dx. Then Pn, n>l, is 

Vi" Jo 

a polynomial of degree n with nonsingular leading coefficient equal to 

1 ^ „ , , ,2,„- 



(1.14) 



r — 2'' 



Moreover, defining Pq 
with respect to W . 



7^ 

1 
2 



7„ = 2+|a|^6"-5n. 
, {Pn)n is a sequence of orthogonal polynomials 



This Rodrigues' formula allows us to find an explicit expression for the polynomials 
{Pn)n in terms of the Hcrmite polynomials. 



Corollary 1.3. For n > 1, we have 

b-'''^-Hn{Vht) ~atb-"/^Hn{Vbt) + fi/„+i(i) 

-2aV'^^nHn^i{Vbt) 2\a\^b"-^^ntHn-i{Vbt) + 2Hnit) J' 

where Hn is the n-th Hermite polynomial defined by Hn{t) — (—1)" (e^* 



(1.15) P„(t) 



In Section 6, using again the Rodrigues' formula ()1.13p . we find the following three 
term recurrence relation for a sequence {^n)n of orthonormal polynomials with 
respect to W. 



Theorem 1.4. The sequence of matrix polynomials defined by = 0, <?3^o 



hVb 







and 



V 71 

V 1 

i^„(i) = A„+i^„+i(t) + Bn^n{t) + Al^n-i{t), n > 0, 



{tt) 4 

(1.16) 

where 

(1.17) 

(1-18) ^ 

are orthonormal with respect to W ll.9\} (where the sequence (7n)n is defined by 

(m i)- 



A„ = y/n 



2bj„ 







\ 


7n-l 




27« 


) 


■f- 


a 


V a 
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This gives for the recurrence coefficients the asymptotic behaviour 

1 

71 



Hm 



= < 




if 6 > 1 



for < 6 < f 



This hmit shows that the recurrence coefficients {An)n do not asymptotically behave 
as a scalar multiple of the identity, as it happens in the examples studied up to now 
in the literature. 



2. Preliminaries 



A weight matrix W is an N x N matrix of measures supported in the real line 
satisfying that (f ) is positive semidefinite for any Borel set A € M, (2) W has 

finite moments / t"dW{t) of every order, and (3) / P{t)dW{t)P* (t) is nonsingular 
if the leading coefficient of the matrix polynomial P is nonsingular (all the matrices 
considered in this paper are square matrices of size N x N). When all the entries 
of the matrix W have a smooth density with respect to the Lebesgue measure, 
we will write W(t) for the matrix whose entries are these densities. Condition (3) 
above is necessary and sufficient to guarantee the existence of a sequence {Pn)n of 
orthogonal matrix polynomials with respect to W. Condition (3) above is fulfilled, 
in particular, when W{t) is positive definite in an interval of the real line. This is 
the case of the weight matrix (|1.8p introduced in this paper. 

The key concept to study orthogonal matrix polynomials (Pn)„ satisfying second 
order differential equations of the form 

(2.1) Pn{t)F2{t) + K{t)Fl{t) + Pn{t)Fo{t) = KnPn{t), 

where the differential coefficients F2, Fi and Fq are matrix polynomials (which do 
not depend on n) of degrees less than or equal to 2, 1 and 0, respectively, is that 
of the symmetry of a differential operator with respect to a weight matrix. Indeed, 
if we write 

i^^(|)V.(*)+(|)^i(t)+Fo(.), 

then D is symmetric with respect to W if and only if the orthonormal matrix poly- 
nomials {Pn)n with respect to W satisfy (|2.ip . where A„, n > 0, are Hermitian 
matrices (see IDlj ). If are not Hermitian, then the operator D can be decom- 
posed as D = Di + iD2, where Di and D2 are second order differential operators 
of the form (jl.ip symmetric with respect to W (see |GT| ). 

The symmetry of a second order differential operator as (|l.ip with respect to a 
weight matrix can be guaranteed by a set of differential equations (which is the 
matrix analogous to the Pearson equation (/2W)' = fiw of the scalar case): 
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Theorem 2.1. For a weight matrix dW — W{t)dt, t E {a,b) (a and b finite or 
infinite) and matrix polynomials F2, Fi, Fq of degrees not larger than 2, 1 and 0, the 
symmetry of the second order differential operator U.l\) with respect to W follows 
from the set of equations 

(2.2) F2W = WF^, 

(2.3) 2{F2WY - FiW ^WF^, {F2W)" - {FiW)' + FqW ^ W F^ , 
under the boundary conditions 

(2.4) 

lim t''F2{t)W{t) ^0, lim t"[{F2{t)W{t)y - Fi{t)W{t)] ^ 0, n>0. 

t~¥a+ ,b- t—^a+,b- 

(See [DGI] or jGPTlj l. 

To check that the weight matrix W (jl.Sp and the coefficients defined by (|1.10p . 
(11.111) and (|1.12p satisfy the differential equations (|2.3|) we will use the following 
theorem: 

Theorem 2.2. Let fl be an open set of the real line. Let F2, Fi, F and T be twice 
differentiable N x N matrix functions on fl, (with T{tQ) nonsingular for certain 
to e fl), and define W{t) — T{t)T*{t). Under the assumptions 

(2.5) F2W^WF^, 

(2.6) T'{t) = F{t)T{t), and 

(2.7) Fi=F2F + FF2+F^, 
we have 

(1) The weight matrix W satisfies the first order differential equation 

2{F2Wy ^ FiW + WF^. 

(2) For a given matrix Fq, the second order differential equation 

{F2W)" - {Fiwy + F„w = wf;, 

holds if and only if the matrix function 

(2.8) x = T-\-FF2F-F'F2~FF^+Fo)T, 
is Hermitian at each point of il. 

Theorem 12.21 is a particular case of Theorem 2.3 of |D3] (the special case when F2 
is a scalar function is Theorem 4.1 of jPGlj ). 



We will need the following theorem (see |D4j ) to find the Rodrigues' formula dis- 
played in Theorem 11.21 for the orthogonal polynomials with respect to the weight 
matrix W (fTOl) : 

Theorem 2.3. Let F2, Fi and Fq be matrix polynomials of degrees not larger than 
2, 1, and 0, respectively. Let W , R„ be N x N matrix functions twice and n times 
differentiable, respectively, in an open set Q, of the real line. Assume that W(t) 
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is nonsingular for t £ and that satisfies the identity S2.2\) . and the differential 
equations i8.3\) . Define the functions P„, n > 1, by 

(2.9) Pn = Rll'^W-\ 

If for a matrix A„, the function Rn satisfies 

(2.10) 



(RnF^)" ~ {Rn [F* + n{F*y]y + i?„ 



F^ + n{F*y+{'l]{F;r 



— A„i?„ 



then the function P„ satisfies 

(2.11) Kit)F2{t) + K{t)Fi{t) + Pn{t)Fa{t) = A„P„(i). 

We will also make use of the following well known formula: for any matrices AT, y G 

inNxN . 



(2.12) ^""^ - {YL^'^xY 



> n>0 



where we use the standard notation 

ad^f Y = Y, &AxY ^ [X, Y]=XY - YX, &d\ Y ^ [X, [X, Y]] , 
and in general, adj^^ Y = [X, [ad^ Y]]. 

3. Symmetric differential operator 

In this Section we prove Theorem ll.il that is, the second order differential operator 
with coefficients given by (|1.10p . (|l.lll) and (|1.12l) is symmetric with respect to the 
weight matrix W (jl.Sp . 

We now list some technical relations which we will need in the proof of Theorem 
11.11 (they will be proved later) . 

Lemma 3.1. Let the function F2 and the matrices A, £/ , 4', & and ^ be defined 
by flJ^) . (T^, (L^, (Eg) and (T^, respectively. Then 

[N]_l 

(3.1) [<^]= E i2j + l)a,A^^+\ 

(3.2) e^** = Fae^*. 

(3.3) e^*^e-^** = -^I - i^_il^e^*^e-^*[^, 

(3.4) *e^*i^e-^* = -^I ~ ^'^^^i^, ^le^^^e"^*. 

r iv-1 1 

(3.5) ^ ^] = t^^r^^A^^- 

^ ' L '^J Z^^ ^(2j + l)J-i 

(3.6) K, - ^ = 2&(jV-\) -^' t-^' ^1 ■ 
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We are now ready to prove Theorem 11.1 



Proof, (of Theorem [H]) 



The symmetry of the second order differential operator with respect to the weight 
matrix W wiU be a consequence of Theorems 12.11 and 12.21 We have to check the 
boundary conditions (|2.4I) and the three equations (|2.2p and 

To make the proof easier to foUow, we proceed in four steps. 



First step: Boundary conditions {2.Ji\j . Proof: Since is nilpotent, we deduce that 
e^^* is a polynomial. The matrix function T = e'^'e®* then decays exponentially 
at oo because the entries of the diagonal matrix S> are negative. Hence the weight 
matrix W — TT* also decays exponentially at oo. Since F2 and -Fi are polyno- 
mials, it follows straightforwardly that rFaM^ and i"[(F2W^)' - FiW], n>0, have 
vanishing limits at 00. 



Second step: F2W = WF^. Proof: Formula of Lemma lO shows that F2T = 
T^P, where 5* is the diagonal matrix (real entries) defined by (|1.5|) and T — e'^'e®* . 
Then 

F2TT* — T'^T* = T{T'^)* = T{F2T)* = TT*F2- 
Since W = TT* , we get that F2W ^ WF*. 

Third step: 2(^2^7)' = FiW + WF^. Proof: To check the equation 2(^2^^)' = 
FiW + WFi, we use the first part of Theorem 12.21 with = M. Hence, we have to 
prove that T satisfies T'{t) ~ F{t)T{t), where F is a solution of the matrix equation 

(3.7) Flit) = F2it)Fit) + F{t)F2{t) + F^{t). 

Taking into account that T = e'^'e®* , a direct computation gives that 

(3.8) F{t) = ^ + 2te-^*^e-^*. 

The definition of F2 (fTTOl) and jSH) give 

F2F + FF2 = + + y -g^K, /] + 2t(«'e^*^e--^* + e^*^e-^**) 

Using the definition of * ([T31) . ([XS]) and ([Q)) of Lemma [Owe get 



F2F + FF2 = 2^ + .Li) (^^ + ^J) + ^] _ 2hlt. 

Formula p.7p now follows easily taking into account the definitions of F2 (jl.lOp . 
Fi (fTTT]) and * (fTS]) . 



foMrt/i step; {F2W)" - (FiVK)' + FqW^ = WF^. Proof: 
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Using 2 of Theorem l2.2l this is equivalent to prove that the matrix 

(3.9) x = T-\-FF2F~{FF2)' + Fo)T 
is Hermitian. 

We actually will prove that the matrix function x defined by (j3.9l) is diagonal with 
real entries. 

Taking into account that T{t) = e'^'e®* , and ^ is diagonal, it is enough to prove 
that the matrix function 

(3.10) e = e-^\~FF2F ~ [FF^]' + Fo)e'^' 
is diagonal. 

We first compute e-'^\FF2F)e'^K 

Taking into account the expression for F{t) in p.Sp . that and e'^* commute and 
using p. 21) . one has after straightforward computations 

The definition of & ([TTB)) and * ([T3)) now give 

(3.11) e-^\FF2F)e'^^ = ^ jr\^ ^ ~ ~ 

We now compute e~-^*(Fi^2)'e-^*. Using again the definition of F ([XS)) and 
of Lemma 13. 1[ we have that 



The definition of 9 (HH) and F2 (fTTO)) give 

(Fi^2)' = (^^* + /] - btl 

Identity p.ip of Lemma \3A] shows that e'^* and [^Z, ^\ commute (they are linear 
combinations of power of A) . One then obtains 

(3.12) e-^\FF2)' e'^' = ^^^K, J] - bL 

We finally compute e--^*fbe-^*. The definition of Fq (fLT2|l gives 

e-^*Foe^* = e-^*(26 J + .s/^-^)e'^' 

= 26e-^*^e^* + j^^e-^'^e^K 

Again p.ip of Lemma l3.1l shows that ^ and ^] commute. Hence ad^ ^ — 0, 
n > 2. Using this fact and (|2.12p one obtains 



The definition of ^ (jl.Sp gives 

(3.13) e"^*Foe^* = 26^ - 2bt[j^ , + ^^^^ / - 
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We now substitute ([3TT|) . (|3?T2l) and ([XTS]) in the definition of £, (fXTOl) obtaining 

p.6p of lemma 13.11 finally gives 

which it is indeed a diagonal matrix. 
This finishes the proof of Theorem 1.1. 

□ 



It remains to prove Lemma 13.1 

Proof, (of Lemma 13. ip 
First step. Proof of i3.1\) : 



Using induction on k, we easily find that [A'', ^\ — kA, k > 1. This shows that 
ad'Xj^ = 0, n > 2. The definition of .o/ ((TTl) gives now (I5T|) . 



Second step. Proof of i2j and {5^: 

First of all, the definition of ^I' (II. 5p gives 
(3.14) K,*] ^""^ 



The definition of £/ (fTTT)) and ((XT|) show that the matrices £/, [£/, J\ and e-^* 
commute (they are linear combination of powers of A). We then have 

(3.15) ad^(^')=0, n>2, 

(3.16) K/]e^* = e^*K,^]. 
Using (Pl^ and (IXTC)) we get 



g.«/t^g-^t = ad". * = + tK, 

n=0 

The definition of * and F2 ([TTU)) give now 
In a similar way, we have 

Using now the definition of '3 (jl.6p . (13.141) and p.l6p we have 

^ 6 (6-a)t 
2 N -\ 

This proves p.3p . The proof of p.4p is similar. 
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Third step. Proof of L3.5\) : Using the definition of ^ in (|1.7p and we have to 

prove the following identity 

i=o i=o j=i ^ ' ^ 

This is equivalent to prove 

[f]-i [t^] 
(3.17) o,„.(2. + l)A^<..>-«. = H^f^ i: o^^jJ^A^. 

j,s=0 j=l w -r ; 

Taking into account that A is nilpotent, 



(3.18) asCtj — aj 



(2(j + .) + l)^+-i 
and writing k = j + s, we find 
(3.19) 



V Q'fe 42(fc+i) V ^ ^ 

^ (2fc+l)'^-l I TO i 

fc=0 ^ ' m=0 ^ ^ 



(2m + l)™(2(fc- m) + 1) 



k — m— 1 



We now use Abel's binomial identity (see for instance [Ro[ p. 18]): for z^w ^ C, 
(3.20) H ( ^ ) ^"^ + ^^"^^ " + w)''^""' = zi;-i(z + w + A:)^ 

m=0 ^ ^ 

Then 

^ ( M (2m + l)'"(2(fc - to) + 1)'=-™-! = 2^(1 + fc)^ 

m=0 ^ ^ 

p.igp now gives 

E «,a.(2.+i)A2(^---)= 5: -.tI^a^'^--^'- 

Using p. 181) one has 

a,+i(j + l)(2j + 1)^-1 1-6 
ttj-ai = — — — , with ai — 



i2j + 3y ' Ab{N-l)' 



Thus, 



'(2j + l)J-" 1-6 f^^ (2j+3)J 
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This proves p.l7p and then p.5|) as weh. 



Four step. Proof of iS.6]) : Taking into account (|3.1I) and p.5|) . we have to prove 
the fohowing identity 

This is equivalent to prove 

(3-21) E 2.o,A^^-+^= E E "^-^T^Trfe^^''^'^''^'- 

Using p. 181) once again, we have 

s=0 j = l ^ 

Writing k — s + j, and taking into account that A is nilpotent of order N, we get 
for the right hand side of p.2ip the expression 

fc=i j=i ^ ' 

Writing m ^ k — j , one obtains 

(2fc + 1)^-1 I ^ ; ^ ^ ^ 2' 



m—l 



k=l ^ ' m=0 

Abel's binomial identity ([3^ now gives 

fc-i 

m / ' 2' 2 



E ( ^ ) - ™)'=-™(m + = {{2k + l)'^ - (2fc + 1)'^-!) 



m— 

From where one can easily deduce p.2ip . 
This proves ([SJ]) . 

The proof of the Lemma is now complete. 

□ 
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4. RODRIGUES FORMULA 



In this section we will prove Theorem 11.21 which provides a Rodrigues' formula for 
a sequence of orthogonal polynomials with respect to the weight matrix W for size 
2x2 (fTOl). 



Write Rn for the functions 
(4.1) i?„ = (-l)"e-*' 



2t 



fe-^ed-")* +^(n + 2t2) 



Erf (\/6t) - Erf (i) 



where, as usual, Erf denotes the error function Eii{z) = —= 



at 
2 

2 

dx. 



The Rodrigues' formula ()1.13|) can then be written as P„ = ^. 

First of all, we explain how one can use Theorem 12.31 to find these functions 
n > 1. 

Indeed, Theorem 11.11 for size 2x2 gives for the weight matrix W the following 
symmetric second order differential operator 

(4.2) 



1 a{b-l)t 
b 



-2bt 




2ab 
-2bt 




2b 



Since the operator D is symmetric with respect to W , the n-th monic orthogonal 
polynomial Pn with respect to W satisfies the differential equation D(P„) = A„P„, 
where the eigenvalues A„ are given by 

-2bn 
-26(n-l) 



A« = 



Theorem 12.31 associates the following second order differential equation (|2.10p to 

the differential operator D (|4.2p : 

(4.3) 



Rn 



1 







a{b-l)t b 



Rn 



-2bt 

a[b{2 + n) - n] 





-2bt 



Rn^n — Rn- 



Take now a solution i?„ of this differential equation and write ¥„ = R^^^W ^. 



Theorem 12.31 guarantees that the function Yn satisfies the differential equation 
D{Yn) = A„y„. Notice that F„ and P„ satisfy the same differential equation. 
We have hence looked for a solution i?„ of the differential equation (j4.3p such that 
the matrix function R^^W"^ is also a matrix polynomial of degree n with nonsin- 
gular leading coefficient. This is the procedure we have used to find the functions 
Rn given by (gH]). 

We now prove Theorem ll.21 which establishes that actually the functions R'h'^W'^ 
define a sequence of orthogonal polynomials with respect to W. 



Proof, (of Theorem [Hg) 
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(n) 

Using the Rodrigues' formula for Hermite polynomials, = (—1)" ( e^* ) e* , 

[Sz| Chapter 5] , we have that 



(-'■) 
(-"■) 

(Erf(t))(" 
(Erf(V5t))(" 



(-1)" 
4 

(-1)" 



{-inVbrHniVbtje 



(-l)"-i6"/2Ai?n-i(V6i)e-''*' 



These identities give, after straightforward computations using the three term re- 
currence relation for the Hermite polynomials <_ff„ = _ff„+i/2 + nHn-i (see |Szl 
Chapter 5]): 



Taking into account that 



e"' -ae"' t \ 



we finally find the expression (jl.lSp for P„ given in Corollarv ll.3l 



-2ab''^^nHn-iiVbt) 2\a\H''^^ntHn-iiVbt) + 2i7„(t) 



This shows that i?i"'(i)W^"i is a polynomial of degree n (note that the entry (1,2) 
of this matrix is actually a polynomial of degree n — 1) with nonsingular leading 
coefficient equal to r„ (|1.14p . 

The orthogonality of P„ and t'^I, < k < n — 1, with respect to W follows taking 
into account that 

' PnW{t)t''dt = J Rl^\t)t''dt, 
and performing a careful integration by parts. 



□ 



5. Three term recurrence relation 



In order to find the recurrence coefficients (jl.l7l) in the three term recurrence rela- 
tion of Theorem 11.41 we have followed the strategy of |DG2j or [PL) . 
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Proof, (of Theorem [L4 



We first compute the norm of the monic orthogonal polynomials P„ with respect 
to W. Using the Rodrigues' formula (|1.13p . we have 

\\Pnf = J Pn{t)W{t)rdt = J Ri:'h"dt, 

where r„ are the leading coefficient of P„ (|1.14p . An integration by parts and the 
formulas for _R„ and r„ (see Theorem 1 1.2p then give 



(5.1) ||P„ 
If we write 

(5.2) A„ = 
the polynomials 



7«+i 
26^ 



V 



ln+1 





2 ' / 



are then orthonormal with respect to W . 

We now prove that the they satisfy the three term recurrence relation (|1.16p . 
This is just a matter of computation. Indeed, the coefficient An+i in (jl.l7p is then 

The formula (|5.2p for A„ gives now the formula for A„ in (|1.17p . 

On the other hand we have for the recurrence coefficient i?„ in (|1.16p the expression 
Bn = A„i?„A,7\ where 

Bn = coeff. of in P„- coeff. of i" in P„+i- 

From (ll.lSp . we get that 



(5.3) Bn = {h + {h- l)n) 

\ 7„7„+i 

and the formula for P„ in (jl.l7p follows easily. 




□ 



The three term recurrence relation for the polynomials (Pn)n (|1.13p can now easily 
be computed. Indeed, taking into account the expression for the leading coefficient 
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r„ (fTM]) of Pn and A„ of we find that P„ = G„^„ where 





(5.4) 



7n+l 



2&- 




/2^ 



In particular, this gives Pq = 



1 
2 



If we write 



(5.5) tP^{t) - i„+iP„+i(t) + B„P„(i) + C'„P„_i(i), n > 0, 

it follows from (|1.16p that A„+i = G„A„+iG^j:j, P„ = G„P„G^^ and G„ = 
G„A*G,7l]^. An easy computation gives now 



1 







An-- 



G„ = n I Irfn 

\ 1 



Q 7n^ , P„ = (-n+(n + l)6) 

7n / 

7«+i 



/ 

2a6"-5 



267n 



7r! + l 



The P^ norm of the polynomials P„ follows easily from the formula (I5.4p for the 
matrices G„: 

/ 7«+i 

|lP„f = 2"0Fn! 2&"+5 

V 27„ 







In a similar way, the three term recurrence relation for the monic orthogonal poly- 
nomials (P„)n can be deduced: 

tP^{t) = Pn+l{t) + BnPnit) + G„P„-l(i), " > 0, 

where P„ is the one in (|5.3p and 

/ 1^ 



V 



In 
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